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Abstract 

The unsteady blade pressures due to the operar 
tion of the propfan at an angle to the direction of the 
mean Sow are obtained by solving the unsteady three- 
dimensional Euler equations. The configuration con- 
sidered is the eight-bladed SR7L propfan operating at 
takeoff conditions and the inflow angles considered are 
6.3° , 8.3°, and 11.3V The predicted blade pressure 
waveforms are compared with in-flight measurements. 
At the inboard radial station ( r/R = 0.68) the phase of 
the predicted waveforms shows reasonable agreement 
with the measurements while the amplitudes are over- 
predicted in the leading edge region of the blade. At 
the outboard radial station ( r/R = 0.95), the predicted 
amplitudes of the waveforms on the pressure surface 
are in good agreement with flight data for all inflow 
angles. The measured (installed propfan) waveforms 
show a relative phase lag compared to the computed 
(propfan alone j waveforms. The phase lag depends on 
the axial location of the transducer and the surface of 
the blade. On ;he suction surface, in addition to the 
relative phase iag, the measurements show distortion 
(widening and steepening) of the waveforms. The ex- 
tent of distortion increases with increase in inflow angle. 
This distortion seems to be due to viscous separation 
effects which depend on the azimuthal location of the 
blade and the axial location of the transducer. 

Introduction 

The enhanced fuel efficiency characteristics of the 
propfans over the conventional propellers were demon- 
strated by NASA through its Advanced Turboprop 
Project. 1 Further improvement in fuel efficiency, noise 
reduction and structures may result only from a better 
understanding of the flow features at the design and 
off-design operating conditions. This understanding is 
crucial to the improvement of design methodologies for 
future propfan design. 

Flight measurements on a large-scale 9-ft (2.74m) 
diameter SR7L propfan have provided a unique 
database of blade stresses, near-field sound pressure 
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levels and blade surface pressures over a range of oper- 
ating conditions. 2 * 3 This database may be used to un- 
derstand the propfan aerodynamic characteristics, val- 
idate prediction methods and improve design method- 
ologies. In the flight tests, the propfan was installed on 
the left wing of a modified instrumented GuHstream II 
testbed aircraft. A nacelle tilt arrangement was used to 
vary the inflow angle to the propfan. Three nacelle tilt 
angles were considered, -3* and -1* (tilt down) and +2* 
(tilt up). The effective inflow angle, a, to the propfan 
however, depends on the airplane angle of attack, na- 
celle tilt and upwash angle of the flow into the propfan. 

The detailed in-flight unsteady blade pressure 
measurements were made on a specially designed in- 
strumented blade using 30 transducers. 3 The blade suc- 
tion surface had 20 pressure transducers distributed 
over three radial stations {r/R — 0.68, 0.86 and 0.95, 
where r is the radial distance and R is the blade tip ra- 
dius) while the pressure surface had 10 pressure trans- 
ducers distributed over two radial stations (r/R = 0.68 
and 0.95). The axial positions of the transducers from 
the blade leading edge are listed in Table I in terms of 
normalized chord length, x/c where x is the axial dis- 
tance and c is the blade chord. Tests were carried out 
over a range of takeoff and cruise operating conditions. 

The unsteady flow features of a propfan at cruise 
operating conditions were examined by N alla sa m y and 
Groeneweg 4 They showed that at inflow angles of the 
order of 5° strong passage shocks extending from suc- 
tion to pressure surface form and dissolve during each 
revolution of the blade. They also studied the unsteady 
flow effects on the blade loading and pressure wave- 
forms at takeoff conditions for an inflow angle of 8.3* 
(corresponding to a nacelle tilt of -1°) 5 , and compared 
with in-flight measurements. The predicted blade 
pressure waveforms showed reasonable agreement with 
flight data at the inboard radial station ( r/R = 0.68) 
whereas at the outboard radial station (r/R = 0.95) 
the agreement was poor. 

In the present paper the effect of inflow angle 
on the blade pressure at takeoff operating conditions, 
Mach number, M = 0.31 and advance ratio, J = 
1.6, is further studied by obtaining unsteady thre*- 
dimensional Euler solution for 6.3° and 113* inflow 



angles corresponding to nacelle tilt angles of -3° and 
+2° respectively, of the flight test. 3 The results ob- 
tained here are analyzed together with those obtained 
for 8. 3°, 5 to understand the effect of inflow angle on 
the blade loading and blade pressure waveforms. The 
predicted blade pressure waveforms are compared with 
flight measurements. 

Numerical Solution of Unsteady, 
Three-Dimensional Euler Equations 

The unsteady three dimensional Euler equations 
governing the in viscid flow through a propfan are solved 
employing a solution procedure developed by Whitfield 
et al. 6,7 In this procedure the Euler equations in conser- 
vative differential form are transformed from a Carte- 
sian reference frame to a body fitted curvilinear refer- 
ence frame. Then the transformed equations are dis- 
cretized employing a finite volume technique. An ap- 
proximate Riemann solver is used for block interface 
flux definitions and a Lower-Upper (LU) implicit nu- 
merical scheme is used to solve the discretized equa- 
tions. The flowfield is represented by a multiblock com- 
posite grid to limit the core memory requirements. 

Flow Configuration and 
Computational Grid 

The configuration considered here is that of the 
eight-bladed SRTL propfan of the flight test. 3 The di- 
rection of rotation of the propfan and the axes of refer- 
ence are shown in Fig. 1. The rotation of the propfan is 
clockwise, looking downstream and the azimuth angle, 
is measured from the z axis as shown. The grid em- 
ployed here to represent the flowfield is the same as that 
shown in Fig. 2 of Ref. 5. An H-gria, with 107 x 41 x 13 
(axial by radial by circumferential) points in each blade 
passage is employed. Each blade passage is divided into 
three blocks with 107 x 41 x 5 grid points in each block. 
Thus 24 blocks of grid were used to describe the entire 
flowfield with 456,248 nodal points. Each blade sur- 
face is represented by 49 x 27 (axial by spanwise) grid 
points with higher resolution near the leading and trail- 
ing edges, the hub and the tip. Similar grids used with 
the present solution technique have produced results 
which are in reasonably good agreement with data.*’ 9 

Results and Discussion 

The unsteady three dimensional Euler solutions 
were obtained for the following two flight test cases of 
Ref. 3: (i) run = 187, id = 3a and nacelle tilt = -3° 
and (ii) run = 113, id = 3a and nacelle tilt = +2°. The 
takeoff operating conditions, Mach number = 0.31 and 


advance ratio = 1.6 are considered. The aircraft angle 
of attack for these test cases was 5.4°. In addition to 
nacelle tilt angle and aircraft angle of attack, the up- 
wash angle at the propfan is needed to determine the 
effective inflow angle. 

In the absence of a simple computational proce- 
dure to determine the up wash angle, an experimental 
correlation obtained by Heidelberg and Woodward 10 
from their SR7A model propfan test in the NASA Lewis 
9 by 15-ft wind tunnel was used to estimate the inflow 
angle to the propfan. They first measured the pres- 
sure response of the blades as a function of inflow angle 
for the propfan alone configuration. Then a wing was 
installed downstream of the propfan (as in the flight 
test case, tractor configuration) and blade pressure re- 
sponse was measured over a range of wing angles of 
attack and nacelle tilt. It is assumed that the local in- 
flow angle and the propfan angle of attack are the same 
in the propfan alone case. Then by matching the mea- 
sured first harmonic response of the blades with wing 
installation to that of propfan alone configuration, they 
found the equivalent inflow angle as a function of wing 
angle of attack (Fig. 13 in Ref. 8). This correlation was 
obtained for takeoff conditions, Mach number s 0.2, 
advance ratio = 0.88, and with a straight wing. The 
same correlation is used here to obtain the effective in- 
flow angle to the propfan. For an aircraft (wing) angle 
of attack of 5.4°, a nacelle tilt (droop) angle of -3° gives 
a value of 6.3° for the effective inflow angle, a, (Fig. 1), 
while a nacelle tilt of +2® gives a value of 11.3®. 

The Euler solutions reported here were obtained 
with inflow angles of 6.3® and 11.3° giving an inflow 
variation of 5° for a 5° change in nacelle tilt. These 
solutions are examined together with that obtained for 
an inflow angle of 8.3° (nacelle tilt = -1°) in Ref. 5, 
at the same takeoff operating conditions. The present 
solutions were obtained from an impulse start for three 
complete revolutions of the propfan, to obtain a rea- 
sonably accurate solution. By the third revolution, the 
results have stabilized as indicated by the periodic vari- 
ation of per blade power coefficient during the second 
and third revolutions of the propfan, Fig. 2. The fig- 
ure shows the variation of the power coefficient with 
azimuth angle for three inflow angles, for a blade start- 
ing at ♦ = 0 and executing three complete revolutions 
(* = 0 to 1080°). The results of the third revolution 
are analyzed and presented here. The predicted total 
power coefficient (for eight blades) are nearly the same 
for all three inflow angles and are about 1% higher than 
the measured value. The expected sinusoidal variation 
of the blade loading due to inflow angle is clearly ob- 
served in all cases. The amplitude of the stabilized 
power coefficient during the third revolution varies +44 
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and -42 percent about the mean for o = 6.3° , +61 and 
-53 percent for a = 8.3°, and +89 and -68 percent for 
a = 11,3°. For the low inflow angle (a = 6.3°) the 
variation about the mean in the positive and negative 
directions is nearly the same, as it was also observed 
in Refs. 5, 9. For higher inflow angles, the difference 
between the minimum and maximum levels about the 
mean results directly from the higher positive inflow 
angles, as a look at the velocity triangles would show. 

The harmonic content of the blade loading is ob- 
tained by Fourier transforming the blade power coeffi- 
cient variation during the third revolution of the blade. 
The power coefficient C p = a 0 + a,* cos art + 6j sin art, 
where a 9 , a,-, 6* are the Fourier coefficients and t is the 
time. The loading spectra are shown in Fig. 3. It is 
seen that the predicted mean blade power coefficient is 
nearly the same for all three inflow angles and the first 
harmonic dominates the blade loading. The first har- 
monic loading lags the blade motion by 3.62°, 3.68* and 
3.75° for the three inflow angles, 6.3*, 8.3* and 11.3*, 
respectively. The relatively low value of lag and the 
smaller variation with inflow angle are due to the low 
reduced frequency and flow unsteadiness at the takeoff 
operating conditions considered here. 

The azimuthal variation of the power coefficient, 
dC p /d{r/R) t are shown in Fig. 4. Figure 4a shows the 
computed elemental power coefficient variation for an 
inflow angle of 11.3° at four azimuthal locations, $ = 
0®, 90* , 180° and 270°. Also shown in the figure is the 
curve for the steady flow (a = 0*). The shape of the 
curve for any azimuthal position is similar to that of 
the steady flow. The magnitude of the cyclic variation 
depends on the spanwise location. The elemental power 
coefficient variation for an inflow angle of 8.3° is shown 
in Fig. 4b. The trends of the curves with azimuth an- 
gle are similar to that of a = 11.3*; higher loading 
levels occur for * = 180* and 270* and lower loading 
levels for $ = 0* and 90* positions. Figure 4c shows 
the elemental power coefficient variations for a = 6.3*. 
The magnitude of the cyclic variation is small for this 
inflow angle. At any inflow angle, the maximum devia- 
tion (in absolute value) from the steady value occurs at 
the radial station where the elemental power coefficient 
attains a maximum, as would be expected from angle 
of attack variations along the blade span. 

The comparisons of the predicted blade pressure 
waveforms with flight data are presented next, for two 
inflow angles, a = 6.3* and 11.3*. Such comparisons 
for a = 8.3° were reported in Ref. 5 and reference to 
it is made where appropriate. In these waveform com- 
parisons, 0* corresponds to the (top-dead-center) ver- 
tical direction for aircraft installation as in the presen- 
tation of flight data in Ref. 3. At the outset, it should 


be mentioned that for the takeoff operating conditk>ns 
considered here, the flight data do not indicate the for- 
mation of a leading edge vortex. During the flight test 
no steady measurements were made. However, some 
steady pressure data were obtained during the unsteady 
pressure test by retaining the DC component of the 
pressure signal. These data do not indicate the for- 
mation of a leading edge vortex. When a leading edge 
vortex is present 9 , the pressure waveforms on the suc- 
tion surface show double hump form which is absent 
in the measured waveforms for the present cases, thus 
suggesting the absence of the vortex(see the waveforms 
an d discussion below). 

Figure 5 shows the unsteady blade surface pres- 
sure as a function of azimuth angle for the transducer 
locations on the suction surface at r/R = 0.68 for 
a — 6.3*. Hie measured and predicted waveform in- 
dicate that the response is the largest near the leading 
edge and the response reduces gradually towards the 
trailing edge. The predicted phases of the waveforms 
are in close agreement with data. However, the ampli- 
tudes are overpredicted, the maximum overprediction 
in absolute value occuring near the leading edge. The 
pressure waveforms for the four transducer locations on 
the pressure surface at this radial station are shown in 
Fig. 6. Here again, the phases of the predicted wave- 
forms agree quite well while the amplitudes are over- 
predicted. 

The pressure waveforms at the outboard radial sta- 
tion r/R = 0.95 for a = 6.3* are shown in Fig. 7 and 
8. Figure 7 shows the pressure waveforms on the suc- 
tion surface. Even at this radial station, the prediction 
shows that the largest response occurs near the leading 
edge and reduce gradually towards the trailing edge. 
The measurements, however, indicate that the maxi- 
mum response occurs at x/c = 0.42. For this and the 
transducer at x/c — 0.58, the magnitudes are severely 
underpredicted. The measured waveform is distorted 
from the linear sinusoidal form of the predicted one. 
The waveform distortion may occur due to viscous ef- 
fects manifested in the form of a small separation bub- 
ble as suggested by the low pressure (within the bub- 
ble) at 4 ~ 180* for x/c = 0.42 and 0.58 and steep rise 
in pressure (outside the bubble) at about ♦ = 270* . 
Hie phases of the waveform also differ from the flight 
data. For all the transducer locations, the measured 
waveform shows a relative phase lag compared to the 
Euler solution. This phase lag also attains a maxi- 
mum at x/c = 0.42. This phase lag seems to be the 
result of an install ation effect, that is, the presence of 
the w ing in the flight tests as compared to the prop- 
fan alone configuration of the computation. Heidelberg 
and Woodward 10 noted similar phase variations in their 
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model tests in the wind tunnel with and without wing 
installation. They show a relative phase lag of 11* with 
the wing installation for a = 10° when compared to 
the propfan alone case at a transducer location r/R = 

0.75 and x/c = 0.1 on the suction surface of the model, 
while the magnitudes are nearly the same. The phase 
lag observed here ranges from 4° near the leading edge 
(x/c = 0.05) to about 60° at x/c = 0.42. On the pres- 
sure side (Fig. 8), however, the magnitudes are well 
predicted, but the relative phase lag exists at ail the 
transducer locations and it ranges from 15° to 50*. It 
should be noted here that the comparison of the wave- 
forms predicted using the same solution technique with 
the propfan alone wind tunnel measured waveforms in 
Refs. 8,9 did not show any such phase lags. 

The pressure waveforms for an inflow angle of 8.3° 
presented in Ref. 5 showed similar behavior at the out- 
board radial station. There the phase lag varied from 
8° to 30° on the suction surface and from 20° to 38° on 
the pressure surface. 

Next, the pressure waveforms for a = 11.3® are 
presented for the two radial stations. At the inboard ra- 
dial station, both the predicted and measured pressure 
waveforms on the suction surface show that the maxi- 
mum response occurs near the leading edge ( x/c = 0.05) 
and reduces gradually towards the trailing edge. At x/c 
= 0.05 both magnitude and phase are in good agree- 
ment with flight measurements. At other transducer lo- 
cations. the amplitude is overpredicted, the maximum 
overprediction occuring at x/c = 0.15. On the pressure 
surface (Fig. 10), the predicted waveforms agree well 
in phase compared to measurements, while the ampli- 
tudes are overpredicted. 

At the outboard radial station r/R = 0.95, the 
pressure waveforms on the suction surface axe shown in 
Fig. 11 for a = 11.3° . First, it is observed that the mea- 
sured waveforms are distorted from sinusoidal form by 
broadening and steepening. This distortion increases 
with increase in inflow angle, 6.3° to 11.3®. At a = 
11.3°, the maximum distortion appears at x/c = 0.25 
transducer location. As noted earlier, the distortion 
of the waveform may be the result of viscous separa- 
tion. The measured waveform at x/c = 0.05 shows a 
relative phase lead compared to the computed (propfan 
alone) waveform, but at all other transducer locations 
a phase lag ranging from 15° to 55® is observed. On 
the pressure surface at this radial station the predicted 
waveforms (Fig. 12) show very good agreement in mag- 
nitude with measurements (except at x/c = 0.92). But 
a relative phase lag of the measured waveforms ranging 
from 7® to 30® is observed. At x/c = 0.92 the blade 
response is low and the predicted waveform is 180® out 
of phase with that of flight measurement. This appears 


to be due to the strong interaction of the blade-wake 
and tip-vortex flows and it is not clear as to the kind of 
interaction that would produce such a phase variation. 

Concluding Remarks 

The unsteady blade surface pressures due to the 
operation of the propfan at an angle to the mean 
flow direction were obtained by solving the three- 
dimensional Euler equations. Three positive inflow an- 
gles, a = 6.3®, 8.3® and 11.3® were considered. The 
predicted waveforms at the inboard radial station show 
reasonable agreement with in-flight measurements. At 
the outboard radial station, the measured waveforms 
show a relative phase lag compared to the measured 
ones. This phase lag is due to the installation effects, 
that is, the installed propfan of the flight test compared 
to the propfan alone configuration of the computation. 
On the suction surface the measured waveforms also 
show distortion (widening and steepening) which in- 
creases with inflow angle. This distortion appears to 
be due to viscous effects which are not considered here. 
Higher flow resolution near the tip region and consider- 
ation of the viscous effects may improve the prediction 
at the outboard stations. 
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Table I Unsteady Pressure Blade 
Transducer Locations in Flight Test 3 


Radius 

r/R 

Surface 

Normalized chord 
r/c 

0.68 

Suction 

.05 .10 .15 .25 .40 .60 .80 


Pressure 

.15 .40 .60 .80 

0.86 

Suction 

.07 .15 .25 .37 .50 .75 .90 

0.95 

Suction 

.08 .25 .42 .58 .75 .92 


Pressure 

.08 .25 .42 .58 .75 .92 
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I'owcr Coefficient, Cp I’owcr Coefficient, Cp 






Fifl. 6 Pressure Waveforms on the Pressure Surface at r/R 0.68, ct 6.3° 
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Rg. 7 Pressure Waveforms on the Suction Surface at r/R=0.95, a-6.3° 


X/c-.OS X/C-.25 




5 


X/c=.42 


X/C«.5fl 



• 5 1 
5| 


X/C*.75 


90 180 270 

Azimuth Angle, Deg. 


seo 



X/C-.92 

Euler Computations 


1 

0 


90 180 270 SCO 

Azimuth AngTe, Deg. 


Fig. 8 Pressure Waveforms on the Pressure Surface at r/R*0.95, a=6.3° 
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Pressure, kPa . Pressure, kPa 





Fig. 9 Pressure Waveforms on the Suction Surface at r/R=0.68, a- 1 1 .3° 
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